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ABSTRACT 

Context. We report on a successful, simultaneous observation and modeling of the sub-millimeter to near-infrared flare emission of 
the Sgr A* counterpart associated with the super-massive (4xlO*Mo ) black hole at the Galactic center. 
Aims. We study and model the physical processes giving rise to the variable emission of Sgr A*. 

Methods. Our non-relativistic modeling is based on simultaneous observations that have been carried out on 03 June, 2008. We 
used the NACO adaptive optics (AO) instrument at the European Southern Observatory's Very Large Telescope and the LABOCA 
bolometer at the Atacama Pathfinder Experiment (APEX). We emphasize the importance of a multi- wavelength simultaneous fitting 
as a tool for imposing adequate constraints on the flare modeling. 

Results. The observations reveal strong flare activity in the 0.87 mm (345 GHz) sub-mm domain and in the 3.8yu/2.2;um NIR. 
Inspection and modeling of the light curves show that the sub-mm follows the NIR emission with a delay of 1.5+0.5 hours. We 
explain the flare emission delay by an adiabatic expansion of the source components. The derived physical quantities that describe 
the flare emission give a source component expansion speed of v„p ~ 0.005c, source sizes around one Schwarzschild radius with flux 
densities of a few Janskys, and spectral indices of q'=0.8 to 1.8, corresponding to particle spectral indices ~2.6 to 4.6. At the start of 
the flare the spectra of these components peak at frequencies of a few THz. 

Conclusions. These parameters suggest that the adiabatically expanding source components either have a bulk motion greater than 
V(,v,, or the expanding material contributes to a corona or disk, confined to the immediate surroundings of Sgr A*. 

Key words, black hole physics, infrared: general, accretion, accretion disks. Galaxy: center, nucleus. Black Holes: individual: SgrA* 



1. Introduction 

At the center of the MiUcy Way, at a distance of only about 8 kpc, 
stellar orbits have convincingly proven the existence of a super- 
massive black hole (SMBH) of mass ~4x10^Mq at the posi- 
tion of the compact radio, infrared, and X-ray source Sagittarius 
A* (Sgr A*; Eckart & Genzel 1996, 1997, Eckart et al. 2002, 
Schodel et al. 2002, Eisenhauer et al. 2003, Ghez et al. 2000, 
2005a). Additional strong evidence of an SMBH at the position 
of Sgr A* comes from the observation of rapid flare activity both 
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in the X-ray and NIR wavelength domains (BaganofF et al. 2001 ; 
Genzel et al. 2003; Ghez et al. 2004, Eckart et al. 2006abc). 
Due to its proximity, Sgr A* provides us with a unique oppor- 
tunity to understand the physics and possibly the evolution of 
SMBHs at the nuclei of galaxies. Variability at radio through 
sub-millimeter wavelengths has been studied extensively, show- 
ing that variations occur on timescales from hours to years (e.g. 
Mauerhan et al. 2005, Eckart et al. 2006a, Yusef-Zadeh et al. 
2008, MaiTone et al. 2008). Several flares have provided evi- 
dence of decaying millimeter and sub-millimeter emission fol- 
lowing simultaneous NIR/X-ray flares. Although much effort 
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Table 1. Observation log. 
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Table 2. Details of NIR observations. 



has been invested in monitoring of the galactic center variability 
at these wavelengths, only a handful of simultaneous observa- 
tions aimed at analyzing their coiTelation have been carried out 
until now (see for instance Eckart et al. 2006a, 2008ab; Yusef- 
Zadeh et al. 2006ab, 2007, 2008, Marrone et al. 2008). 

Here we report on the first measurements of Sgr A* that 
successfully detected simultaneous flare emission in the near- 
infrared and sub-millimeter domain using the ESO VLT and the 
APEX sub-mm telescope^]. In the second section we describe 
the observations and data reduction followed by an outline of 
the flare modeling in section 3. We conclude with a discussion 
and summary in sections 4 and 5. 

For optically thin synchrotron emission we refer throughout 
this paper to photon spectral indices (a) using the convention 
Sy oc V " and to spectral indices (p) of electron power-law dis- 
tributions using N{E) oc with p - (I + 2a). The assumed 
distance to Sgr A* is 8 kpc (Reid 1993), consistent with recent 
results (e.g. Lu et al. 2008, Ghez et al. 2005a, 2008, Eisenhauer 
et al. 2003). 



2. Observations and data reduction 

The sub-mm regime is of special interest for simultaneous flare 
measurements, as it provides important constraints on the mod- 
eling. In this general wavelength range synchrotron source com- 
ponents that also radiate in the infrared domain become optically 
thick, and represent the dominant reservoir of photons that are 
then scattered to the X-ray domain through the inverse Compton 
process. On 3 June 2008, substantial progress was made in this 
matter: as part of a larger campaign, Sgr A* was simultaneously 
observed in the sub-millimeter and NIR wavelength domains us- 
ing the European Southern Observatory facilities. These obser- 
vations result in the first simultaneous successful detection of 
strongly variable emission in both wavelength domains. In the 
following subsections we describe the data acquisition and re- 
duction for the individual telescopes. The basic details of the 
individual observing sessions are given in Table[T] 



' Based on observations with the ESO telescopes at the Paranal 
Observatory under programs IDs:077.B-0028, 79.B-0084, and 81.B- 
0648. The sub-mm data are based on observations with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fiir Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory 



2.1. The NIR data 

Near-infrared (NIR) observations of the Galactic center (GC) 
were carried out with the NIR camera CONICA and the adaptive 
optics (AO) module NAOS (briefly "NACO") at the ESO VLT 
unit telescope 4 ( YEPUN) on Paranal, Chile, during the night be- 
tween 2 June and 3 June 2008. In all observations, the infrared 
wavefront sensor of NAOS was used to lock the AO loop on the 
NIR bright (K-band magnitude ~6.5) supergiant IRS 7, located 
about 5.6" north of Sgr A*. Details on integration times and ap- 
proximate seeing during the observations are listed in Table|2] 
Here A is the central wavelength of the broad-band filter used. 
DIT is the detector integration time in seconds. NDIT is the 
number of exposures of integration time DIT that were aver- 
aged on-line by the instrument. N is the number of images taken. 
The total integration time amounts to DITxNDITxN. Seeing is 
the value measured by the Differential Image Motion Monitor 
(DIMM) on Paranal at visible wavelengths. It provides a rough 
estimate of atmospheric conditions during the observations. The 
atmospheric conditions (and consequently the AO correction) 
were stable during the observations. 

All observations in the L'-band (3.8 fim) and K-band 
(2.2 fim) were dithered to cover a larger area of the GC by mo- 
saic imaging. The sky background for the K-band observations 
was extracted from the median of stacks of dithered exposures 
of a dark cloud - a region practically empty of stars - a few ar- 
cminutes to the northwest of Sgr A*. The L'-band observations 
were interspersed with frequent off-target exposures that served 
to determine the rapidly varying background in this filter. All ex- 
posures were sky subtracted, flat-fielded, and corrected for dead 
or bad pixels. Subsequently, PSFs were extracted from these im- 
ages with StarFinder (Diolaiti et al. 2000). The images were 
deconvolved with the Lucy-Richardson (LR) algorithm. Beam 
restoration was caiTied out with a Gaussian beam of FWHM cor- 
responding to the respective wavelength. The final resolution at 
2.2 and 3.8 ;um is about 60 and 104 milli-arcseconds (mas), re- 
spectively. 

The flux densities of the sources were measured by aperture 
photometry with circular apertures of 52 mas radius and cor- 
rected for extinction, using Ak - 2.8, and Au - 1.8. The rel- 
ative flux density calibration was caiTied out using known K- 
and L'-band flux densities of IRS16SW, IRS16C, IRS16NE, and 
IRS21 by Blum, Sellgren & Depoy (1996). This results in K- and 
L'-band flux densities of the high velocity star S2 of 22+1 mJy 
and 9+1 mJy, respectively, which compare well with magnitudes 
and fluxes for S2 quoted by Ghez et al. (2005b) and Genzel et al. 
(2003). The relative photometry for Sgr A* was done with the 
known fluxes and positions of 9 sources within V!6 of Sgr A*. 
The measurement uncertainties for Sgr A* were obtained on the 
reference star S2. The background flux in the immediate vicin- 
ity of Sgr A* was obtained by averaging the measurements at 
six random locations in a field located about 0'.'6 west of Sgr A* 
that is free of obvious stellar sources. 
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(mJy) 
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(mJy) 


I 


NIR L'-band 


<04:43 


06:20+10 mm. 


05:15+5 min. 


>90 


40+5 


32+2 


>29 


<3 


II 


NIRL'-band 


06:20±10min. 


07:20±10min. 


06:50+5 min. 


60+10 


40+5 


17+2 


>14 


<3 


III 1 


NIR L'-band 


07:20+10 min. 


>8:20 


08:10+5 min. 


>55 


>25 


17+2 


>14 


<3 


III 2 


NIR Ks-band 


<8:20 


08:54+10 min. 


08:43+5 min. 


>20 


10+5 


12+1 


>6 


<6 


IV 


NIR Ks-band 


08:54+10 min. 


09:48+10 min. 


09:20+5 min. 


55+10 


22+5 


14+1 


> 8 


<6 




sub-mm 


<06:00 


08:50+20 min. 


07:00+15 min. 


>220 


>140 


4250+200 


850+200 


-3400 



Table 3. Here we list data of the flares observed by the VLT and APEX on 3 June 2008. We list the estimated start and stop 
times, the full width at zero power (FWZP) and full width at half maximum (FHWM) values, as well as estimates of the flare peaks 
and light curve minima which may be considered as offsets. The flare peak flux densities are given in 10"^ Jansky. Flux density 
uncertainties are given as their Icr values. 
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Fig. 1. The infrared light curve of Sagittarius A* on 3 June 2008. The 
data are represented by vertical red bars (+lo") and a black connection 
line between them. Upper panel: L'-band flare events. For a better rep- 
resentation we added a constant value of 50 mJy to the light curve for 
the reference star Lower panel: The K-band flare emission immediately 
following the L'-band observations. In this case no offset was added to 
the light curve for the reference star. 



The resulting light curves of Sgr A* and S2 (as a constant 
reference star) are shown in Fig.[T] 



2.2. The sub-mm data 

The Galactic center 870 jum data were taken with the newly 
commissioned LABOCA bolometer array, located on the 
Atacama Pathfinder Experiment (APEX) telescope at the Llano 
Chajnantor, Chile, at an altitude of 5105 m. The observations 
were carried out as part of global multi-wavelength monitoring 
program in May/June 2008. The radiation is collected by the 
APEX telescope and fed through conical horns to the detectors 
of the Large APEX Bolometer Camera (LABOCA), an array of 
295 composite bolometers that are extremely sensitive to con- 
tinuum radiation. With a total bandwidth of about 60 GHz the 
system is optimized for the 345 GHz atmospheric window. 

The full width half maximum (FWHM) of the point spread 
function (PSF) at 870 /im is ~ 19". To allow for the best possible 
reconstruction of the complex sub-mm emission in the GC re- 
gion we used on the fly (OTF) maps perpendicular to the galactic 
plane (Weiss et al. 2008, Siringo et al. 2007). The maps were ob- 
served with a scanning speed of 37s modulating the source sig- 
nals even for extended structures into the useful post-detection 
frequency band of LABOCA (0.1-12.5 Hz). To avoid scanning 
artefacts in the reconstruction of the extended thermal emission 
surrounding Sgr A*, we derived and averaged maps from the 
data with an inclination of -10°, 0° and H-10° with respect to an 
axis orthogonal to the galactic plane. The mapping steps per- 
pendicular to the scanning direction was 30" and the integration 
time 280 seconds, yielding an rms noise level of ~ 150 mJy/beam 
for each map. 

Immediately after each galactic center map, either G 10.62 
or IRAS 16293-2422 were observed as secondary calibrator 
sources. We reduced the data with the BoAQ software package. 
The data reduction process included correction for atmospheric 
zenith opacity (towards the Galactic center t~0.3 up to 0.7), 
flat-fielding, de-spiking, correlated sky noise removal, and the 
removal of further correlated noise due to instrumental effects. 

After correcting for pointing offsets between individual maps 
(determined from the position of Sgr A* in each map) all scans 
were co-added to obtain high S/N, 48'x25' fully sampled maps 
of the sub-mm emission in the GC regions. From this combined 
map the point source Sgr A* was modeled with a Gaussian and 
taken out. The uncertainty in that process will result in a small 
offset of the sub-mm light curve. We estimate that this effect is 
less than 10% of the peak flux density of SgrA*. In this way 
we created a model image of the extended 870//m emission sur- 
rounding Sgr A* with Sgr A* effectively set in an off state (see 
Fig.©. 



^ BoA:'http://www.astro.uni-bonn.de/boawiki/Boa] 
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Fig. 2. Left: Sagittarius A* at the beginning and the peak of the APEX measurements. Right: A section of the larger map of the Galactic center. 
The location of Sagittarius A* is indicated by a circle. Sagittarius A* has been subtracted out as this map is used as a reference for the extended 
structure. The colors code the flux density in Janskys. 



Each data point of the Sgr A* sub-mm light curve (see 
Fig.lSj was derived from the model subtracted maps, modeling 
a Gaussian source and deriving the peak. In order to provide a 
stable and consistent calibration we used the t values derived 
from the sky dip and radiometer measurements, as well as sec- 
ondary flux calibrators and a comparison to the reference map. 
Between the individual observing epochs the flux densities of 
the secondary flux calibrators vary by <10%. From a compari- 
son of different reference sources, we estimate the relative point- 
to-point uncertainty in the calibration of the light curve measure- 
ments to be of the order of 4% (1 cr). The remaining flux density 
variations in the light curve in Fig.|3]may arise from residual fast 
and uncompensated opacity or pointing changes. 



3. Flare Analysis 

3.1. Adiabatically expanding source components 

Our basic assumption to model the sub-mm light curves is the 
presence of an expanding uniform blob of relativistic electrons 
with an energy spectrum n(E) oc E^'' threaded by a magnetic 
field. As a consequence of the blob expansion, the magnetic field 
inside the blob declines as R^^, the energy of relativistic particles 
as ' and the density of particles as R^^ (van der Laan 1966). 
The synchrotron optical depth at frequency v then scales as 



^ y(p+4)/2 I ^ x-(2p+3) 



TO - 



Rc 



and the flux density scales as 



Roj 1 - exp(-To) 



(1) 



(2) 



Here Rq, So, and To are the size, flux density, and optical depth 
at the peak frequency of the synchrotron spectrum vo. The goal 



of the present model is to combine the description of an adia- 
batically expanding cloud with a synchrotron self-Compton for- 
malism, as this is the most likely physical scenario to explain the 
delay between the sub-mm and the simultaneous near-IR and X- 
ray peaks. Thus, we use the definition of tq as the optical depth 
coiTesponding to the frequency at which the flux density is a 
maximum (van der Laan 1966), rather than the definition of to 
as the optical depth at which the flux density for any particular 
frequency peaks (Yusef-Zadeh et al. 2006b). This implies that to 
depends only on p through the condition 



-to(;'-h4)/5- 1 =0 



(3) 



and ranges from, e.g., to 0.65 as p ranges from 1 to 3. 
Therefore, given the particle energy spectral index p and the 
peak flux 5o in the light curve at some frequency vo, this model 
predicts the variation in flux density at any other frequency as a 
function of the expansion factor (R/Rq). 

Finally, a model for R{t) is required to convert the depen- 
dence on radius to time: we adopt a simple linear expansion at 
constant expansion speed Ve^p, so that R - Rq -^exp {t - k))- For 
t < to we have made the assumption that the source has an op- 
tical depth that equals its frequency dependent initial value Ty 
at R = R(). So in the optically thin part of the source spectrum 
the flux initially increases with the source size at a constant Ty 
and then decreases due to the decreasing optical depth as a con- 
sequence of the expansion. For the ~4xlO^M0 super-massive 
black hole at the position of Sgr A*, one Schwarzschild radius 
is Rj=2GM/c^~10'" m and the velocity of light corresponds to 
about 100 Rj per hour. For f > to the decaying flank of the 
curve can be shifted towards later times by first, increasing the 
turnover frequency vo or the initial source size Ro, and second, 
by lowering the spectral index asy„ch or the peak flux density So- 
Increasing the adiabatic expansion velocity v^xp shifts the peak 
of the light curve to earlier times. Adiabatic expansion will also 
result in a slower decay rate and a longer flare timescale at lower 
frequencies. 
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Fig. 3. The sub-mm light curve of Sagittarius A* on 3 June 2008. The 
data are represented by vertical red bars {±\a in length) and a black 
connection line between them. 



3.2. Modeling the flare emission 

The combined L' - and K-band data show violently variable emis- 
sion with at least 4 prominent flare events (I-IV, see Fig.lTJ. As 
we switched from L- to K-band during an ongoing flare we at- 
tribute the L-band flare section III.l and the K-band flare sec- 
tion III. 2 to the same NIR flare III with a duration from about 
7:30 UT to about 9:00 UT. 

As can be seen from the first data points in Fig. [3] the sub- 
mm data start with a low flux density at the time of the first 
NIR flare (Fig.[IJ. The gap shown in the sub-mm data between 
5 h and 6 h occurred during culmination of Sgr A*, when the 
source rises above the elevation limit (80°) for observations with 
APEX. In the sub-mm domain the variable emission is dom- 
inated by a single flare that is wider than the individual NIR 
flares (see Table [3j. The so far six reported coordinated SgrA* 
measurements that include sub-mm data (Eckart et al. 2006b; 
Yusef-Zadeh et al. 2006b; Marrone et al. 2008, including this 
work) have shown that the observed submillimeter flares follow 
the largest event observed at the shorter wavelengths (NIR/X- 
ray; see detailed discussion in MaiTone et al. 2008). We therefore 
assume that the sub-millimeter flare presented here is related to 
the observed IR flare events. 

Separating the variable sub-mm emission: In order to sep- 
arate the rapid intra-day flare emission from the sub-mm emis- 
sion that varies much slower on time scales of days, a constant 
flux density level has been subtracted. In most instances Sgr A* 
has been tracked over long time spans on each day to obtain 
light curves. We have determined the median and median devi- 
ation of the lowest flux density measurements of the obtained 
light curves. This allows us to determine a flux density contribu- 
tion that is constant or only slowly variable on time scales much 
longer than a single day. For this purpose we compared 5 APEX 
light curves taken with LABOCA between May 27 and June 
3 (centered at a wavelength of 0.87 mm, Garcia-Marin 2008, 
in prep.) as well as 6 CSO measurements at 0.85 mm wave- 
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Fig. 4. The infrared (top) and sub-mm (bottom) light curve of 
Sagittarius A*. The data are represented by vertical red bars i±icr) 
and a black connection line between them. The dashed line represents 
a smoothed version of the data (after application of a 7 point sliding 
average for all data points exept the first). The models as described in 
the text and in Tab.|4](blue, black, red for models A, B, C, respectively) 
are shown as solid lines (blue: ~2.5 THz; black: ~2 THz; red: ~1 THz 
cutofi" frequency in the source component spectra.) 



lengfli by Yusef-Zadeh et al. (2006a, 2008) and 2 SMA light 
curves at 1.3 mm wavelength by Marrone et al. (2008). The 
SMA interferometer resolves out extended emission on scales 
of a few arcesonds. In the case of the APEX measurements a 
model of the extended emission on scales of 19" is subtracted. 
The flux density offsets of the different light curves are very sim- 
ilar. Therefore we attribute them to be intrinsic to the emission 
of Sgr A*. The comparison results in a constant offset value of 
3.0+0.5 Jy. For the the lowest submillimeter flux densities ob- 
tained on June 3 we find a constant contribution of 3.4 Jy. This 
value is in very good agreement with the median derived above. 
Therefore we subtracted it from the light curve in order to extract 
the intra-day variable part of the sub-mm data. 

Adiabatic expansion modeling: We model the 4 flare events 
I-IV and the smoothed version of the sub-millimeter data 
with 4 adiabatically expanding source components (see dashed 
line in Fig.|4] The calculations show that the initial synchrotron 
self-absorption cutoff frequency of the source components is a 
critical parameter that allows us to distinguish between sepa- 
rate models. From our modeling we can constrain the cutoff 
frequency between 1 and 2.5 THz (model A:2.5 THz, model 
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B:2 THz and model C: 1 THz), and establish from inspection 
and modeling of the light curves a delay between the sub-mm 
and NIR emission of 1.5±0.5 hours. Model fits to the variable 
sub-mm emission are shown in comparison to the data in Fig.|4] 
The corresponding model parameters are given in Tab. |4] To 
select the intra-day variable part of the sub-mm light curve a 
flux density of 3.4 Jy has been subtracted from the sub-mm 
data in Fig.|4] This amount is attributed to more extended (many 
Schwarzschild radii) source components. We attribute the time 
difference between the NIR and sub-mm flares to an adiabatic 
expansion of synchrotron source components with an expansion 
speed of about 0.5% of the speed of light (1500 km/s). While the 
sub-mm light curve is modeled by the evolving optically thick 
part of the spectra, the NIR light curve is modeled through the 
optically thin part of the spectra. We find that the flares are as- 
sociated with source components that have initial sizes of the 
order of one Schwarzschild radius and spectra that peak around 
2 THz. The flux densities are a few Janskys and the spectral in- 
dices are optically thin between the sub-mm and the infrared 
with 0.6<Q'<1.5. Models with significantly different expansion 
speeds or source sizes fail to represent either the extent or shape 
of the observed flare features. 

Models A and C represent border cases: For higher cut- 
off" frequencies than the one used in model A, the required 
flare fluxes at the THz turnover frequencies become increasingly 
higher to account for the sub-mm emission after adiabatic expan- 
sion. At the same time, the spectral index becomes increasingly 
steeper in order to meet the NIR flux densities. Towards cutoff' 
frequencies well above 2.5 THz the spectral index will become 
steeper than what has been measured in the NIR (Eisenhauer et 
al. 2005, Hornstein et al. 2006, Gillessen et al. 2006, Krabbe et 
al. 2006). On the other hand the mean spectral index for the com- 
ponent in model C is closer to the value of 0.6+0.2 obtained by 
Hornstein et al. (2007), but the predicted source component la- 
beled / shows already an unacceptable discrepancy between the 
flux density observed with APEX at 04:56 UT. This discrepancy 
increases with decreasing cutoff frequency. Our favored model 
B is presented in Fig.|5] It has an intermediate cutoff frequency, 
shows an overlap in parameters with models A and C, and con- 
sistently reproduces the flares in the sub-mm and NIR wave- 
length domains. With a deviation of about 3-4cr, the agreement 
of the spectral indicies with results by Hornstein et al. (2007) 
lies between that of model A and C. The spectral indices are 
in good agreement with that obtained for aMmix-mv ~ 1-3 by 
Eckart et al. (2004, 2006a). This is consistent with the fact that 
the optically thin synchrotron spectral index (sub-mm to NIR) 
is expected to equal the broad band spectral index of the SSC 
spectrum. 



SSC modeling as an additional constraint: We use a syn- 
chrotron model with an optically thin spectral index a and rel- 
ativistic electrons with jg ~ 10^ following Gould (1979) and 
Marscher (1983). In the context of Sgr A* this model has been 
described by Eckart et al. (2004, 2008a). This formalism allows 
us to estimate the Synchrotron Self Compton (SSC) X-ray flux 
densities and to compare them with values that have been mea- 
sured before. In addition it allows us to estimate the SSC contri- 
bution of the NIR flux density and the magnetic field. In our 
adiabatic expansion modeling we have assumed that the NIR 
flux density is dominated by the synchrotron emission of the 
THz peak expanding synchrotron component. A strong excess in 
SSC NIR emission would violate that assumption and would re- 
quire modeling the (unfortunately) unknown X-ray flux density 
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Fig. 5. Our favored model B compared to the infrared (top) and sub- 
mm (bottom) light curve of Sagittarius A*. The data are represented by 
vertical red bars (+lcr) and a black connection line between them. For 
furher explanation see caption of Fig.|4] 



as well. The magnetic field can be obtained via B ~ ff^v^^S^^ 
as a direct result of the synchrotron source parameters. Previous 
models of the Sgr A* emission have resulted in magnetic field 
strengths of a few up to 70 Gauss (e.g. Eckart et al. 2006a, 2008a, 
Yusef-Zadeh et al. 2008, Marrone et al. 2008). The source com- 
ponents assumed for the adiabatic expansion model should not 
violate these results. 

Our SSC calculations, that reveal typically observed X-ray 
flux densities in the range of a few 10 to 100 nJy with mag- 
netic fields of a few 10 G, indicate that the best solution lies 
between models B and C in Tab. |4] For model A we find that the 
SSC contribution of the NIR flux density becomes comparable 
to the synchrotron contribution. A correspondingly lower pre- 
dicted NIR synchrotron contribution would require even steeper 
spectral indices. The SSC calculations, however, already indi- 
cate that the spectral indices and cutoff frequencies in model A 
need to be lower by Aq'~0.2 and Av of a few 100 GHz in order 
to give reasonable solutions that agree with typically observed 
and derived X-ray flux densities and magnetic fields. Therefore 
model A clearly appears less preferable. 

Quality of the fit: Sub-structures deviate by up to 4 times the 
measurement uncertainties from the overall structure described 
by flares I to IV. In order to estimate the fit quality we scale 
the assumed uncertainties to within which we can model the 
NIR data to Act. To compare the NIR and APEX data with the 
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adiabatic expansion models, we calculated reduced values. 
Models with a large number of flare or source components would 
increase the number of free parameters and allow us to obtain 
better reduced values for the fits However, more complex 
models are currently difficult to justify. Following the ratio of 
the 41 APEX to 208 NIR data points, we increased the weight 
of the APEX data by a factor of 4. With a total of 4 components 
(each with a time to, peak flux So, size Ro, and spectral index), 
a common cutoff frequency, sub-mm flux offset and expansion 
velocity v^ we obtain a total of 19 free parameters. As shown 
in the Figs.[l]to |5]the L- and K-band data were compared to the 
L- and K-band model predictions, respectively. For our favorite 
model we calculated the uncertainties of the model parameters 
by varying them until the reduced ;^-values increase by A;^'~0.5. 
This corresponds to an increase across either of the 4 individ- 
ual flare intervals in the sub-mm or the NIR domain of A;^~ 1 .0, 
corresponding to a confidence level of about 68.5% (i.e. a 1 cr 
uncertainty) for the model parameters for each of the flare com- 
ponents. 

3.3. Discussing fhe modeling results 

The adiabatic expansion results in a time difference between 
the peaks in the VLT and APEX light curves of about 1.5+0.5 
hours, and compare well with the values obtained in a global, 
multi-wavelength observing campaign by our team in 2007. 
Back then, two bright NIR flares were traced by CARMA 
(Combined Array for Research in mm-wave Astronomy; 
100 GHz) in the US, ATCA (AustraHa Telescope Compact 
Array; 86 GHz) in Australia, and the MAMBO bolometer at the 
IRAM 30m in Spain (230 GHz; first results given by Kunneriath 
et al. 2008, Eckart et al. 2008b). This light curve complements 
our parallel 13, 7, and 3 mm VLBA run (Lu et al. 2008). Other 
recent simultaneous multi-wavelength observations also indicate 
the presence of adiabatically expanding source components 
with a delay between the X-ray and sub-mm flares of about 
100 minutes (Eckart et al. 2006a, Yusef-Zadeh et al. 2008, 
Marrone et al. 2008). As pointed out by Eckait et al. (2008ab), 
a combination of a temporary accretion disk with a short jet can 
explain most of the properties associated with infrared/X-ray 
Sgr A* Hght curves (Eckart et al. 2008ab). 

The adiabatic expansion model: The May 2007 polarimet- 
ric NIR measurements (Eckart et al. 2008a) showed a flare event 
with the highest sub-flare contrast observed until now. In the 
relativistic disk model these data provide evidence of a spot 
expansion and its shearing due to differential rotation. An ex- 
pansion by only 30% will lower the Synchrotron-Self-Compton 
(SSC) X-ray flux significantly. In the framework of the spot 
model this flare event provides additional support for expansion 
of individual source components. 

The expansion speed: The rapid decay of the observed 
NIR/X-ray flai-es (e.g. Baganoff et al. 2001; Genzel et al. 
2003; Ghez et al. 2004, Eckart et al. 2006abc), as well as the 
current results from coordinated observing campaigns including 
sub-mm monitoring (Eckart et al. 2006b; Yusef-Zadeh et 
al. 2006b; Marrone et al. 2008 and this work) suggest that 
non-radiative cooling processes, such as adiabatic expansion, 
are essential, although the adiabatic cooling model results 
in very low expansion speeds. From modeling the mm-radio 
flares Yusef-Zadeh et al. (2008) invoke expansion velocities 
in the range from Vej:p=0.003-0.1c. This compares well with 
the expansion velocity of the order of 0.0045±0.0020 c that 



we obtain in our case with the June 2008 data (see Table |4]i. 
These velocities are low compared to the expected relativistic 
sound speed in orbital velocity in the vicinity of the SMBH. 
The low expansion velocities suggest that the expanding gas can 
not escape from Sgr A* or must have a large bulk motion (see 
discussions in Marrone et al. 2008 and Yusef-Zadeh et al. 2008). 
Therefore the adiabatically expanding source components either 
have a bulk motion larger than M^^p or the expanding material 
contributes to a corona or disk, confined to the immediate 
surroundings of Sgr A*. An expansion of source components 
through shearing due to differential rotation within the accretion 
disk may explain the low expansion velocities. The recent 
theoretical approach of hot spot evolution due to shearing is 
highlighted in Eckart et al. (2008a) and Zamaninasab et al. 
(2008; see also Pechacek et al., 2008) 

Structure of the light curve: The flux density variations of 
Sgr A* can be explained in a disk or jet model (see e.g. discus- 
sion in Eckart et al. 2006ab, 2008a), or they could be seen as a 
consequence of an underlying physical process that can mathe- 
matically be described as red-noise (Do et al. 2008, Meyer et al. 
2008). In our case we find a light curve structure that consists of 
maxima separated by about 70 to 110 minutes with additional 
fluctuations of smaller amplitude. Assuming the presence of a 
disk and by simultaneous fitting of the previously obtained light 
curve fluctuations and the time-variable polarization angle, we 
have shown that the data can be successfully modeled with a 
simple relativistic hot spot/ring model (Meyer et al. 2006ab, 
2007, Eckai-t et al. 2006ab, 2008ab). In this model the broad 
near-infrared flares (~100 minutes duration) of Sgr A* are due 
to a sound wave that travels around the SMBH once. The sub- 
flares, superimposed on the broad flare, are thought to be due 
evolving (i.e. expanding) hot-spots that may be relativistically 
orbiting the central SMBH. The spot emission would then be 
due to transiently heated and accelerated electrons which can 
be modeled as a plasma component (scenarios in which spiral 
wave structures contribute to the observed variability are also 
under discussion, e.g. Karas et al. 2007). 

In case of a jet (see e.g. Markoff, Bower & Falcke 2007, 
Markoff, Nowak & Wilms 2005), the observed flux density vari- 
ations may more likely be a result of the variations in the accre- 
tion process (or jet instabilities) - possibly followed by an adia- 
batic expansion of the jet components - rather than being a result 
of a modulation from an orbiting spot. In this case one may ex- 
pect that red-noise variations on these short times scales are a 
natural extension of the variability found for longer periods. 

4. Summary and discussion 

We have presented new simultaneous measurements of the near- 
IR sub-mm flare emission of the Sgr A* counterpart associated 
with the SMBH at the Galactic center. The data were obtained 
in a global campaign carried out on 3 June 2008, using NACO 
at VLT and LABOCA at APEX. 

The highly variable near-IR emission presents four major 
events (I - IV), whose peaks are separated by about ~80 min. 
The peak of the sub-mm light curve is delayed with respect to 
that of the near-IR data. We argue that this delay is due to the 
adiabatic expansion of synchrotron source components that be- 
come optically thin. The expansion velocity is about 1500 km/s 
(0.5% of the speed of light). We model the four flare events with 
4 source components. In this case the sub-mm light curves of 
the individual components are blended, with a time delay for the 
peaks of about 1.5±0.5 hours with respect to the near-IR data. 
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model flare to 
label label hours 



"exp 

in c 



[Jy] 



(^synch 



Rn 



vo 
[THz] 



reduced 
value 



I 

n 

m 

rv 



5.3+0.2 
6.9+0.5 
8.2±0.4 
9.5+0.2 



0.0040±0.0020 



11.4±3.6 
4.8±3.5 
5.6±3.0 
4.1±2.0 



1.77±0.20 
1.67±0.10 
1.77±0.10 
1.45+0.10 



1.1±0.10 
0.7±0.40 
1.0±0.70 
0.6+0.50 



2.50±0.45 



2.64 



I 

II 
III 

rv 



5.3±0.2 
6.9+0.5 
8.2±0.4 
9.4+0.2 



0.0045±0.0020 



7.8+1.5 
1.4±1.5 
3.3±1.6 
2.8±0.5 



1.57±0.10 
1.25+0.10 
1.56±0.10 
1.33+0.10 



1.3±0.30 
0.9+0.60 
1.2±0.70 
0.9+0.50 



2.00+0.45 



1.82 



I 

n 

m 

rv 



5.3±0.2 
6.9±0.5 
8.2+0.5 
9.4+0.2 



0.0053±0.0020 



4.0±1.3 
1.3+1.0 
0.5±1.3 
1.4±0.3 



1.15+0.30 
1.05±0.20 
0.85+0.20 
1.00+0.20 



1.3+0.60 
1.0±0.50 
1.2+0.30 
0.7+0.40 



1.00+0.30 



2.32 



Table 4. Source component parameters for the adiabatic expansion model of the 3 June, 2008 flare. The flare times to are given 
with respect to the peak of the brighter NIR flares labeled I, II, III, IV. For the adiabatic expansion velocity v^p, the optically thin 
spectral index Osynch and the cutofl" frequency vo the uncertainties were derived over the entire data set. In addition to v^p the Ro 
values are responsible for the position and width of the infrared flares peaks in time. Different values for Us^.nch are requred to match 
the infrared flux densities. The derivation of uncertainties and reduced values is described in the text. Model B is the preferred 
model (see text). 



The light curve structure is consistent with the previously 
found variability. It could be interpreted as emission from ma- 
terial in relativistic orbits around the SMBH. However, variable 
emission from a jet or an explanation as a dominant flux den- 
sity contribution from a red-noise process cannot be excluded as 
well. 

These data show that the VLT/APEX combination is espe- 
ciaUy well suited for very long simultaneous light curves be- 
tween the NIR and the sub-mm domain. Further simultane- 
ous radio/sub-mm data, NIR K- and L-band measurements in 
combination with X-ray observations should lead to a set of 
light curves that will aUow us to prove the proposed model 
and to discriminate between the individual higher and lower en- 
ergy flare events. Here Chandra's high angular resolution is ide- 
ally suited to separate for weak flares the thermal non-variable 
bremsstrahlung and the non-thermal variable part of the Sgr A* 
X-ray flux density. 
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